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We introduce a general formalism to describe the effects of Markovian noise which is spatially 
correlated, typically decaying over some finite correlation length. For any system of interest, this 
formalism describes spatial correlations without the necessity to choose a particular microscopic 
model for the environment. We present a method of mapping the equations to Lindblad form and 
discuss functional forms for homogeneous spatial correlation functions. We also discuss two exam- 
ple microscopic models for the environment which exhibit non-trivial spatial-temporal correlation 
functions. 
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The study of open quantum systems and the concept 
of a density matrix master equation underlies much of 
modern quantum physics, be it quantum optics, atom 
optics, condensed matter physics or quantum computa- 
tion. When deriving or assuming a particular form for 
the master equation via the usual system-bath modelp]], 
it is common to assume each component of a system ei- 
ther couples to the same bath {correlated or collective 
decoherence channels) or individual baths (uncorrelated 
or independent decoherence channels). We explore the 
regime between these two extremes, introducing the con- 
cept of a correlation length £ and deriving a general mas- 
ter equation method for treating such partially correlated 
environments. 

Decoherence induced by a correlated environment is a 
commonly observed effect. The concepts of super- and 
sub-radiance 0, [j| and decoherence free subspacesj3"ll| 
rely on the interference resulting from several subsystems 
coupled to the same environmental mode. Even some of 
the foundational work on decoherence in quantum com- 
putation considered both fully correlated and uncorre- 
lated environments Recently, the ramifications of 
correlated environments have been discussed in such di- 
verse situations as scalable quantum error correction ll3r- 
Il9j |. photosythesis and biological chromophores (20l424| 
and multi-atom trapping experiments 25h 2a | . 

The Lindblad equation^], [2^-31 1 is the workhorse of 
open- quantum systems due to its simple form and well 
behaved mathematical properties. Yet, deriving a gen- 
eral master equation of this form for (partially) correlated 
environments is non-trivial. Lindblad operators can be 
derived or assumed which act individually, pair-wise or 
collectively, yet how does one choose these in a physically 
sensible manner? We use a general Bloch-Redfield ap- 
proach, where environmental noise correlation functions 
appear naturally in the formalism. Given sufficiently well 
behaved environmental correlations, a closed form master 
equation can be obtained with the same form as the Lind- 
blad equation but whose operators and rates are linked 
directly to the original physical system/bath Hamilto- 
nian. Following this route we consider a generalisation 
of the environmental correlation function which includes 



spatial (as well as temporal) correlations and therefore 
derive a general master equation. This describes spa- 
tially correlated decoherence independent of the particu- 
lar bath Hamiltonian and purely based on environmental 
correlation functions. We also consider several examples 
where such spatial-temporal correlation functions can be 
derived microscopically. 

Starting from the usual system-bath Hamiltonian, 



H = H S + H L 



Hi 



(1) 



comprised of system (Hs) 
between them 



bath {Hb) and the interaction 



Hint — 



(2) 



where Sj are system operators and Bj bath operators. An 
arbitrary basis {|a„)} and the Hamiltonian eigenstates 
Hs |wn) = w n \<jJn) are connected via the transformation 
matrix V = Yl n \uj n ) (a n \. A compact and general form 
of the Bloch-Redfield equations is then given by (see sup- 
plementary material): 



p= n [p 



with 



-pVqjkVhj + SjpVqjkV^ (3) 



(a n \qjk\a m ) = (a„\ V ] s k V \a m ) -C jk (oj m - ui n ) (4) 
(a n \q jk \a m ) = (a n \ V 1 ' s k V \a m ) -C kj (uj n - ui m ) (5) 



C jk (uj) 



dre iur (4 (r)B fc (0)) 



(6) 



The spectral functions Cj k (u>) define the bath, i.e. the 
environment. They are given by the Fourier transform 
of the correlation functions of the bath operators in the 
interaction picture Bj(r) = exp(iHBT)Bj exp(— ?'i7eT). 
The spectral functions define both the correlation func- 
tion with increasing time difference r and which of the 
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pairs of bath operators Bj , B k (corresponding to system 
operators Sj,Sk) are correlated. 

The operators Bj can be grouped such that 



Cjk{uj) if j and k are in the same group 
if j and k are in different groups 



This may be due to either the assumption of a certain 
structure in the environment or groups of different types 
of coupling operators (see supplementary material). We 
call each group an independent bath since they are not 
correlated with the other baths. 

If the environment is thought to be a continuum with 
decaying correlations over increasing distances a more 
flexible model is to introduce a spatial dependency in 
the spectral function: 



C jk (u),rj,T k ) := 



dre™ (^(T,r,)B fc (0,r fe )) 



where the components of the system couple to the envi- 
ronment at positions rj, v k . This function is similar to 
a van-Hove- function [32], however we use real-space coor- 
dinates rather than k-space. As in conventional Bloch- 
Redfield theory, this spectral function can be derived 
from a more fundamental microscopic model or its form 
can be phenomenologically assumed based on the physi- 
cal nature of the environment in question. 

Strictly speaking the spectral function Cj k {u) is ac- 
tually given by a one-sided Fourier transform Dj k {u>) — 
Cjfe(w) + iFjk(uj), i.e. may be complex. Such a term 
would lead to additional coherent dynamics. In the sec- 
ular approximation it can be written as a correction 
Hg Hg + H cor to the system Hamiltonian (see supple- 
mentary material), where H cor = Y]j k Fjks\.Sj. In cases 
of uncorrelated decoherence (all terms for which j ^ k are 
zero) this correction of the system Hamiltonian is gener- 
ally neglected since the terms s^Sj in qubit systems are 
diagonal in the system eigenbasis. In cases of correlated 
decoherence however the correction Hamiltonian H cor for 
qubit systems can lead to interaction terms between the 
qubits, i.e. environmentally induced interactions. 

The Bloch-Redfield equations do not guarantee com- 
plete positivity by their mathematical form but depend 
on a consistent model of the environment's temporal 
and spatial correlation functions. When assuming a 
phenomenological spatial-temporal correlation function 
Cj k {ui, Tj, r k ), highly non-trivial physical conditions ap- 
ply for multipartite correlations. We therefore present a 
way to check for a mapping to the Lindblad equations 
which are known to be "the most general form of the 
generator of a quantum dynamical semigroup" [l , 2^, 3(| . 
This means that the time-evolution operator due to a 
certain Markovian master equation is a completely posi- 
tive map of the density matrix if and only if there exists 
a set of operators in which the master equation takes on 
Lindblad form. 



We take the Bloch-Redfield equations in the secular 
approximation (see supplementary material): 



P = Ti {p.H,i 



jk s k ps 



\{ S \^P]) ( 7 ) 



This equation can be mapped to Lindblad form if and 
only if there exists an invertible matrix W such that 
W'CjkW is diagonal with non-negative eigenvalues. The 
eigenvalues then correspond to the Lindblad rates. The 
Lindblad operators are given by: 



For real- valued homogeneous correlation functions which 
depend only on distance, C jk (rj,r k ) = O jk (\rj - r k \), 
the coefficient matrix Cj k is symmetric and its diagonal- 
isation is unitary, i.e. eq. [7| can be mapped to Lindblad 
form if the coefficient matrix is positive semi-definite. If 
one or more eigenvalues are negative then according to 
Sylvester's law of inertia[33| all transforms W^Cj k W will 
have the same number of negative values on the diagonal, 
i.e. there is no mapping to Lindblad form. 

To illustrate this procedure, we consider a one- 
dimensional chain of equidistant subsystems (separated 
in space by a distance d). The correlations in the en- 
vironment felt by these subsystems decay exponentially 
with distance according to Cj k = cxp(— a\j — k\) where 
the parameter a > sets the decay in terms of the envi- 
ronmental correlation length, a = We prove this ex- 
ponential form can always be mapped to Lindblad form. 

We first note that homogeneous coefficient matrices 
Cjk with row j and column k are Toeplitz-matrices. To 
determine whether it is positive semi-definite we utilise 
Lemma 6 in Chapter 4.2 of reference [HI], which applied 
to our notation reads: 

The eigenvalues of an nxn Hermitian Toeplitz matrix 
in the notation 



( t i_i 

ti to 



-n+l 



\ 



(8) 



\ tn-l ■ ■ ■ to J 

with absolutely summable elements t m , 

oo 

y] \t m \ < oo 

m— — oo 

are not smaller than the minimum and not greater than 
the maximum of the Fourier series /(A) defined by 



/(A)= £ t m e 
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For the case of exponentially decaying correlations, 
this Toeplitz matrix is specified by the elements u m = 
e -o|m|_ yj a a geometric series, they are found absolutely 
summable: 



E 



2^e 

ni—1 



1 - e~ a 



1 < oo 



Similarly the Fourier series evaluates to: 

oo 

g(X) = u m cxp(-imX) 

2(1 - e" cosA) 



771 — — OO 



-f 



f - 2e~ a cosA + e 



-2a 



> 



Since the Fourier series 17(A) is greater than zero, the 
eigenvalues of the corresponding Toeplitz matrices are 
greater than zero. Any master equation based on an 
exponentially decaying spatial correlation function Cjk = 
exp(— a\j — fc|) can therefore be mapped to Lindblad form. 

Another useful example is that of Gaussian decay, 
Cjk = cxp[— a(j — fc) 2 ], in which case the elements of 
the corresponding Toeplitz matrix, v m = e _al ™ 2 , are ab- 
solutely summable, 



E ki = 2 E 



1 < 00 



m— — cxd 



m=0 



while the Fourier series is given by, 

00 00 
h(X) = v k e lkX = 1 + 2 e ~ ak2 cos(fcA) > (9) 



fc— — 00 



fc=i 



which can be shown using Jacoby's ima gina ry transfor- 
mation (chapter 21, page 475 of reference |35j). Again we 
see that the Fourier series h(X), and therefore the eigen- 
values of the corresponding Toeplitz matrices are greater 
than zero which means that any master equation based 
on a Gaussian correlation function can be mapped to 
Lindblad form. 

Finally we give a counterexample where the mapping 
fails, indicating an unphysical correlation function. The 
Toeplitz matrix corresponding to a step correlation func- 
tion has the elements: 



iv, , 



1 for |m| < 2 
otherwise 



The elements are absolutely summable X)m=-oo \ Wr - 
3 however the Fourier series has negative values: 

1 = 1 + 2 cos A 



J2 w k e lkx = e lX + e~ iX 




In fact for the specific case of a system consisting of three 
subsystems 



(10) 



i.e. the matrix is not positive semi-definite. This 3x3 ma- 
trix is a principal submatrix[3(| of the coefficient matrix 
for any correlation function which corresponds to a step 
function. Followin g th e Definiteness Criterion of sym- 
metric matrices (36l. |37|. it follows that no step function 
correlation matrix (other than trivial limiting cases) is 
positive semi-definite. Step functions should therefore 
not be used as the form of a homogeneous spatial corre- 
lation function as they always lead to a time evolution 
which is not completely positive, even though the tem- 
poral component can be considered to be Markovian. 

To illustrate the concept of spatially correlated deco- 
herence, we now consider two example microscopic mod- 
els where the environmental space-time correlation func- 
tion can be calculated explicitly. In both examples we 
see how the concept of correlation length arises naturally 
from the interplay between thermal effects and the prop- 
agation of excitations within the environment. 

One of the most studied system-bath models is the 
spin-boson modelf^, [12, 38|-40| where a spin undergoes 
dephasing due to an ensemble of non-interacting bosons. 
As an example of an environment bath with non-trivial 
correlations, we consider a one-dimensional chain of in- 
teracting bosons, the 'tight-binding chain', described by 



N 



N-l 



ata. 



X x+l C 



(11) 



Transfor ming into fc-space via a lattice Fourier 
transform[41( diagonalises the Hamiltonian for large N, 
Hb = QJk a k a k with Uk = uj q — 2g cos kd where d is the 
lattice spacing. For a finite chain length, the sum over k 
is given by k n = 2?:n/Nd for integer n £ [— N/2, N/2]. 

Without choosing a particular system Hamiltonian we 
choose the bath operators Bj in the interaction Hamil- 
tonian (eq. [2]) to be the lowering and raising operators 



at a particular position x, 



ikxd 



a k , 



B\x) 



E 



p—ikxd 



N 



(12) 



fc— — 00 



which are non-local in fc-space [4l|. If we assume the 
steady-state density matrix of the harmonic oscillator 
chain is a fully mixed state (in the Fock basis of the 
k-modes) we can compute the four spatial-temporal cor- 
relation functions for a coupled chain of harmonic oscil- 
lators: 

(B(t,x)B(0,x?)) = (B^{t,x)B^{0,x')) = (13) 
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and 



(£t( T) X )B(0, x')) = £ ^e~^-^ d e^ T (n k ) (14) 

k 

(B(t,x)&(0,x')) =^2—e lk ^ x "> d e- tUkT (l + (n k )) 



(15) 



where (n k ) is the expectation value of the occupation 
number of mode k. These expressions are exact and can 
be evaluated numerically for finite N. To obtain more 
physical insight we take the limits of large energies, where 
(n k ) w exp(— huJkP) with /3 = 1/fc^T, and a long chain 



N/2 

E 

n = -N/2 



" *n n 2lT J -it Id 



(16) 



If we assume that the correlations in the chain are dom- 
inated by propagating excitations, we may also linearise 
the dispersion relation, u k w Lo q + 2</(|fc|d — tt/2) within 
the first Brillouin zone. 

Taking these limits, we find that the correlation func- 
tion decays for large r as 



lim (B\t,x)B(0,x')) oc 



1 



2irgT 



(17) 



The bath therefore causes "Markovian noise" if and only 
if the coupling strength g in the chain is larger than the 
characteristic energies in the system. 

For linearised dispersion, we can also evaluate the spa- 
tial correlations at one moment in time t = t' r = 0, 
in the limit of g(3 3> and normalised by the value for 



(&(t,x)B(t,x')) = 
(Bf(t,x)B(t,x)) (2/?ff) 2 + (x - x'f 



(18) 



This indicates that the spatial correlations decay over a 
characteristic correlation length £ = 2d(3g, due to the 
interplay between excitations hopping along the chain 
and thermal noise. 

The spectral function of eq. [TUcan also be computed 
for this system in the limit of large N, using the iden- 
tity S (f(k)) = J2 3 \f'(k/) \ wnere /(fcj) = 0- For cosine 
dispersion and a Boltzmann distribution for (n(u>)) we 
find 

C B i B (uj,x-x') =0(2|ff| - |w + w,|) 

Tx-x' (^) HH)) 



TT^Ag 2 - (lj + u; q y 



(19) 



where T n (x) are the Chebyshev polynomials of the first 
kind. The spectral function shows that the noise consists 
of frequencies centered around the oscillator's eigenener- 
gies ui q with a noise cut-off at u q ± 2g due to the band gap 
of the chain. The spectral function C BB \ (ui, x — x') cor- 
responding to eq. [TS] takes the same form with u> — > —lu 
and (n(M))-).(l + (n(M))). 

A second microscopic example we consider is when the 
system dephases due to the influence of a classical one 
dimensional Ising chain of N coupled spins S x = ± 1 with 
Hamiltonian: 



JV-l 



Hb — — J SxSx+l 



(20) 



with coupling strength J. The spatial correlations in a 
thermal equilibrium state are given by 42 1 



{S X S X >) = tanhl x - x 'l(^J) = exp {ln[tanh(/3J)] \x - x'\} 

(21) 



in our notation. This can be extended [43|,|44( to a space- 
time correlation function by introducing a switching rate 
per unit time a/2 between the two states S x = ±1 for 
each spin (due to a heat bath): 



(Sx(0)S x >(t)) = J2 V lx - x ' +ll Ii(ia\r\) 



(22) 



/ — — oo 



where r\ = tanh(/3J); 7 = tanh(/?2J) and I n (x) is the 
modified Bessel function of the first kind. Note that set- 
ting r — > in eq. [22] yields back eq. [2T] 

Dephasing in qubit systems is generally caused by the 
noise at zero-frequency. The spectral function (eq. [5]) in 
this environmental model at uj = can be calculated 
(assuming positive J) by integrating over each summand 
and evaluating the sum afterwards, 



C(u = 0,x-x') 



2 (|a: - x'\ + C) Cv lx ~ x ' 1 



(23) 



with £ = cosh(2J/3). Here we see the dephasing rate de- 
cays with distance \x — x'\ where the effective correlation 
length depends on Jf3. 

We have presented an extension of the Bloch-Redfield 
formalism which includes spatial correlations within a 
decohering environment. This approach is very general 
as the effect of the environment is modelled via spatial- 
temporal correlation functions. We also provide a proce- 
dure for mapping to Lindblad form master equations as 
well as testing for physicality of phenomenological cor- 
relation functions in the Markovian regime. Two mi- 
croscopic examples were presented where the spatial- 
temporal correlation function can be derived explicitly. 
In these examples, the characteristic correlation length 
depends on the coupling strengths and temperature of 
the environment. 
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DERIVING MARKOVIAN MASTER EQUATIONS FOR SPATIALLY CORRELATED DECOHERENCE: 

SUPPLEMENTARY MATERIAL 

The derivation of the Bloch-Redfield equations [45- 47 1 is a well known procedure which can be found in several 



good textbooksfH l3ll |48|, In this supplementary material, we summarise this derivation and its associated results 
for two reasons, firstly to clarify our notation and form of the equations and secondly to discuss several important 
approximations explicitly. Although these approximations are included in the usual discussions, the introduction of 
spatial correlations requires careful treatment of these points. 

Derivation of the matrix form of the Bloch-Redfield equations 

From the usual system-plus-bath approach: 

H = H s + H B + Him with H mt = SjBj 

j 

where the interaction is given as a product of system operators Sj and bath operators Bj one finds the standard 



master equation in the Born-Markov-approximation (see e.g. Carmichael[31| eq. 1.34) 



P = J2 f°° dt' {[sj^Sk^p - s k (t')psj(t)] (Bi(t)B k tf)) + IpSkWSj® - ~ Sj (t)ph(t')} (Bk^Bjit))} 
ih J o 

where the tilde implies the interaction picture. To define our matrix notation in an arbitrary basis {|a„)} we introduce 
the transformation matrix V = ^2 n \ui n ) (a n \ to the system Hamiltonian's eigenstates Hs \ui n ) = w« \uj n )- To perform 
the t' integration explicitly, we transform back to the Schrodinger picture. For the first term it reads explicitly 

s 3 V y^e- tHsT s k e lHsT Vy ^p(t) (24) 

with t = t — t'. Inserting unity operators J2 n l^n) gives: 

Q = V \ai) (ui\ e- %HsT \u m ) (u m \ s k \u n ) (co n \ e lHsT \u p ) (a p \ (25) 
, " * ' s * ' 

Irnnp 

= J2 lom) M su |w„) (a„| e i ( w -- w ") T (26) 

mi i 

(a m | Q |a„) = (u m \ s k \u n ) e l ^~^ T = (a m | VU k V \a n ) e*^-—-.^ (27) 

Performing these operations analogously on all four terms we can then integrate element-wise and we find a compact 
and general form of the Bloch-Redfield equations: 

P=jj\p, H ^ + ^2 J2 (sjVq jk Vip + Vq jk V^p Sj - P Vq ]k VU + s jP Vq jk V^ (28) 

with (a n \q jk \a m ) = (a„| V*s k V \a m ) D jk (uj m - cj n ) (29) 
(a n \q jk \a m ) = (a„\ Vh k V \a m ) D* jk {u n - ui m ) (30) 

pOC 

D jk {u)= dr e iu)T (Bj(r)B k (0)) (31) 



If Dj k (u>) contains an imaginary term, this leads to additional coherent dynamics (which in the secular approximation 
can even be written as a correction of the system Hamiltonian (see section "Complex spectral functions"), therefore 

only the real part of Dj k causes decoherence. Assuming the property (Bj(— r)B k (0)) = ({Bj(T)B k (0))\ this real 

part can be rewritten in terms of the Fourier transform: 

Re(D jk ) = \c ]k (32) 
C jk (u)= dr e l " T (B 3 {T)B k (Q)) (33) 
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The Markov approximation, which is used to derive the Bloch-Redfield equations, requires Cj k (ui) to be constant on 
the frequency scale of the quantum system, i.e. "smoothness". This means a fast decay of the temporal correlation 
function, which is independent of the spatial correlation length. 



Secular approximation of the Bloch-Redfield equations 



The secular approximation can be used to simplify the Bloch-Redfield equations to the form of eq. [7J First we sep- 
arate the system operators Sj into several operators by multiplication with projection operators onto the Hamiltonian 
eigenstates \u> n ). 

S j'( e ) = E S 3 l Wm ) ( Wm ' ( 34 ) 

The sum extends over all Hamiltonian eigenvalues uj n and ui m with a fixed difference e. Assuming all Sj to be Hermitian 
we find 



Using Eq. I34[ we write: 



Qjk 
qjk 



Sj (-e)= Sj (ey (35) 

E^( £ ) ( 36 ) 

e 

V^8 k (e)±C jk {e)V (37) 

V*52s k (e)±C kj (-e)V (38) 



Replacing eq. [36] to [38] in the Bloch-Redfield equations (eq. [3]) and changing to the interaction picture: 



e iHst Sj {e)e~ lHst = e~ let S] (e) (39) 

One then finds: 

P = E (sji^s^p + s k (e)p Sj (e')) + C kj (-e) (-ps k (e) Sj (e') + s 3 (e')ps k (e))} 

jkee' 

Since e includes negative frequencies the only non-oscillating terms are for e' = — e. In the secular approximation, 
the oscillating terms which are neglected are fast oscillating on the time scale of the system dynamics. This time 
scale ts = 1/| e — e'| for e ^ e' must be much shorter than the decoherence time scale set by the terms Cj k and the 
magnitudes of the Sj. Roughly this can be restated as H$ ^> Hi nt . Applying the secular approximation one finds: 

p = ^yCj k {e) (-Sj(-e)s k (e)p + s k (e)psj(-e)) + C kj (-e) (-ps k (e)Sj(-e) + s i (-e) / 5s fc (e))] ^ 

jke 

Since the sum over e extends over all positive and negative energy differences we can replace e — > — e in the second 
line of eq. [40j Note that this makes the last two terms the Hermitian conjugate of the first two terms. Furthermore 
we swap the equivalent indices j and k in the last two terms. This together with eq. 1351 yields: 

'P = E E \ C ^) (2a*(e)pa J -(c) t - {sj(e)U k (e),p}) (41) 

e jk 

In several cases one finds only one frequency in the elements of qj k and cjj k in which case the summation over e 
contains only one summand and can be omitted, simplifying eq. [41] (back in the Schrodinger picture) to: 

p = l -[ P , H a ] + ^ e \ c * ( 2fl fcp«5 - i s h^ p}) ( 42 ) 

j,k 
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For ease of notation, the form of Eq. (|4"2"j) is used without loss of generality. For cases where the summation over 
e in eq. [41] is relevant, each frequency component is added independently and all further results apply analogously. 
In practice the secular approximation for a particular system is best applied via the method given in appendix Al of 
reference 49fl- The general arguments above are analogous to the derivation of the Lindblad form on page 128 ff. in 
reference If . 

In general, the coupling between a system and its environment can be subdivided into two components, a transverse 
component (where energy can be exchanged between system and environment) and a longitudinal component (where 
the environment simply imparts a phase shift on the system). In the secular approximation only these two coupling 
types (transversal and longitudinal) need to be considered and the two coupling types can be approximated as coupling 
to independent baths. Therefore typically no mixed couplings of transversal and longitudinal operators appear. 

Complex spectral functions 

Strictly speaking the spectral function Cjk{ui) is actually given by a one-sided Fourier transform (see eq. I3ip 
Djk(ui) = Cjk{u) + iFjk(uj) which can be in general complex. Taking that into account eq. HOI becomes: 



Z 5 = E D Jk( e ) (- s j(- £ ) s fc( e )P + Sk(t)pSj(-6)) + Y^( D jk{e))* (-ps k (-e)sj(e) + Sj(e)ps k (-e)) 

jke jke 

Dividing the expression into its real and imaginary components, one finds additional terms to eq. 1411 



(43) 



P = EE (2* fc (e)p* 3 -(e)t - { Sj (e)U k (e),p})+J2Jl iF ^ K^*^) " ^V^(e)p) (44) 

e jk e jk 

which (back in the Schrodinger picture) can be regarded as a correction to the system Hamiltonian, often referred to 
as Lamb shift terms (lj. 

p = j[p,H s + H cor ) + —J2J2~ C ik(e) (2s fc (e)ps 3 -(e)t - {s 3 f(e)+s fc (e), p}) 



e jk 



E jk 

Again for ease of notation the summation over e is omitted in the main text. 
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